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Abstract. Liquid 3 He at low temperatures is an ideal substance to study because 
of its natural purity. The superfluid state, which appears at temperatures below 3 
mK, has many unusual and exciting properties. These have extensively been studied 
experimentally, and the results are in many cases well understood theoretically using 
quasiclassical approximation. In this work we study the case that liquid 3 He is inside 
of aerogel. Aerogel is a very porous substance, where 98% of the volume can be empty. 
The purpose is to understand how the properties of the superfluid are modified when 
the quasiparticles are scattered from an impurity, such as aerogel. We study extensively 
the homogeneous scattering model. From it we derive the Ginzburg-Landau theory of 
impure superfluid 3 He. We give expressions for measurable quantities and compare 
them with experiments. We consider random anisotropy and its effect on the NMR 
properties. More sophisticated scattering models are briefly discussed. 1 



1 Introduction 

Liquid 3 He at low temperatures is the purest substance in Nature because usual 
impurities simply will fall down by gravity. This is one of the reasons why liquid 
3 He, in spite of its strong particle-particle interactions, is one of the best un- 
derstood condensed matter systems (Vollhardt and Wolfle 1990). Especially the 
superfluid state, which occurs at temperatures below 3 mK, shows many com- 
plicated but still well understood phenomena. But the purity can also be a limi- 
tation because often important physical effects exist only because of impurities. 
For example, except of a few marginal cases, all pure elemental superconductors 
are of type I. That is, type II superconductivity, where magnetic field penetrates 
into the superconductor as quantized flux lines, exist because of impurities which 
radically change the properties of pure elements. 

One type of impurity is formed by the surfaces of the container that hold the 
liquid. There has been several studies of liquid 3 He in different confined geome- 
tries, for example between parallel plates or in packed powders (for example, 
Tholen and Parpia 1992). In these cases the effect on 3 He arises from walls, 
i.e. from 2-dimensional interfaces between 3 He and a foreign object. It is very 
difficult to support zero-dimensional (point-like) impurities in liquid 3 He. How- 
ever, there is the intermediate case of one-dimensional impurities. This type of 
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Fig. 1. The measured superfluid transition points of 3 He in aerogel: Porto and Parpia 
(1995, A), Sprague et al. (1995, ■), and Matsumoto et al. (1997, o). The dashed lines 
denote the normal-superfluid and the A-B transition lines in pure 3 He. The A and B 
phases also refer to pure 3 He. 



impurity can approximately be realized by silica aerogels. These are very porous 
materials where, for example, 98% of the volume is empty. Several recent ex- 
periments investigate 3 He in aerogel (Porto and Parpia 1995, Matsumoto et al. 
1997, Sprague et al. 1995, 1996). It is found that, similar to the bulk liquid, 
3 He in aerogel goes to the superfluid state, but both the superfluid transition 
temperature T c and the amplitude of the superfluid state are reduced compared 
to the pure case. The measured suppression of T c is shown in Fig. 1. 

In this article we present some theoretical ideas concerning 3 He in aerogel. 
We first describe the structure of aerogel (Section 2). The main effect of aerogel 
on superfluid 3 He arises from scattering of the quasiparticles of 3 He from the 
aerogel impurity. This effect can be calculated to a good approximation using 
the quasiclassical theory. The general assumptions of quasiclassical scattering 
models are discussed in Section 3. In Section 4 we study a model where the 
scattering is considered as homogeneously distributed. We will limit to tem- 
peratures in the neighborhood of the transition temperature T c , and derive the 
Ginzburg-Landau theory. We calculate measurable quantities and compare them 
with experiments. We find that more sophisticated scattering models are needed 
in order to explain the measurements quantitatively (Section 5). The anisotropy 
of the scattering is considered in Section 6. There we study a random field model 
using similar arguments as presented by Imry and Ma (1975), and find that the 
NMR properties crucially depend on the anisotropy. 
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2 Aerogel 

A nice introduction to aerogels is given by Fricke (1988). The fabrication uses a 
gelation process in the liquid phase. A variety of aerogel materials are possible, 
but the present ones are of silica, (SiC^n- In order to get aerogel it is crucial 
to preserve the gel structure when the solvent is removed. This is not trivial 
because the surface tension of the liquid-gas interface would make the gel to 
collapse in straightforward drying. Fortunately, it is possible to go continuously 
from the liquid to the gas phase. For that one has to move along such a path 
that goes around the critical point (To, po) in the temperature-pressure plane. 
In order to avoid too high temperatures and pressures, methanol (512 K, 8 MPa) 
is preferred to water (647 K, 22 MPa) as a solvent. 

The experiments with superfluid 3 He use samples where the aerogel fills only 
2% of the total volume, V = 0.02. The surface to volume ratio is determined 
by measurements with 4 He, and they give A — 260000 cm -1 (Kim ct al. 1993). 
Based on these numbers alone, we can make some estimations of the dimen- 
sions in aerogel. If we assume that the structure is a lattice of one-dimensional 
strands, we calculate the strand diameter 4V/A = 3 nm and the distance be- 
tween neighboring strands V '4ttV /A = 20 nm. The mean free path for straight 
line trajectories is estimated as i = 4/ A = 150 nm. Although this picture is cer- 
tainly very idealized, the numbers are not very different from the ones obtained 
by other means such as transmission electron microscopy or low-angle x-ray 
scattering. The distance between strands certainly has large variation because 
of fluctuations in the local aerogel density. 

The next step in understanding superfluid 3 He in aerogel is to study the 
characteristic lengths of 3 He. The smallest scale is the average distance between 
neighboring 3 He atoms. This is on the same order of magnitude as the Fermi 
wave length Af = 2w/kf « 0.7 nm. This is somewhat smaller than the diameter 
of the strands. A much larger scale is the coherence length £ = fivf/ 2n kQT c0 . 
Here T c0 is the superfluid transition temperature in bulk 3 He and Vf is the Fermi 
velocity. As a function of pressure, £o decreases monotonically from 77 nm (zero 
pressure) to 16 nm (3.4 MPa, the solidification pressure). The coherence length 
is on the same order as the average strand spacing, but taking into account 
fluctuations in strand density, we can well expect voids of size larger than £o- 
Finally, the mean free path of quasiparticles in pure 3 He is 1 ^m or more at the 
temperatures we are interested in (T < 3 mK). Thus the scattering from strands 
is the dominant factor that limits the mean free path in aerogel. 

3 Quasiclassical Scattering Models 

Before going into the details of different scattering models, it is useful to study 
the general assumptions. Large part of the ideas presented here are given in 
a more mathematical form by Buchholtz and Rainer (1979). The following is 
applicable to both metals and 3 He at low temperatures. Here a "particle" refers 
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to an electron in the former case and to a 3 He atom in the latter. "Superfluidity" 
denotes both the superconductivity of metals and the superfluidity of 3 He. The 
many-body hamiltonian of the system can be written as 



Here H puic contains the kinetic energies of the particles and interactions between 
them, and also interaction with the periodic crystal lattice in the case of metals. 
It may also contain the interaction with external magnetic field, for example. The 
latter term is the impurity potential. In the case of metals, a common source for 
it is that some atoms of the lattice are of different type than the others. In the 
case of 3 He, such a term is introduced by the aerogel. Quite generally, the total 
impurity potential V can be written as a sum so that each term Vj is nonzero 
only in some small region of space. These regions are referred to as "scattering 
centers" . 

The general quasiclassical theory for H purc is discussed extensively by Serene 
and Raincr (1983), and also elsewhere in this book. The validity of this theory 
for the superffuid state relies on the assumption Af <C £o> which is well satisfied 
in 3 He and in most superconductors. The result is that instead of strongly inter- 
acting particles, the system is better understood in terms of quasiparticles. The 
quasiparticles propagate through the medium similar to classical particles. The 
interactions between the quasiparticles are weak, and can be neglected in many 
cases, which would not be the case for particles. Such quasiclassical description 
is valid for properties that are dominated by processes taking place near the 
Fermi surface in the momentum space. This includes virtually all phenomena in 
the superfluid state. 

Let us consider adding the impurity effects to the quasiclassical theory. We 
will not set a limitation to the magnitude of the impurity potential i.e. it can 
be large, on the order of the Fermi energy Ef. The crucial assumption is that 
this large potential occurs only in a limited volume of the total space. This 
means that the quasiparticle description is still valid in most of the space. The 
average distance a quasiparticle can travel between it collides with scattering 
centers defines a mean free path I. In order to use the quasiparticle picture, 
£ has to be large in comparison to the Fermi wave length Af. In addition, we 
neglect coherent scattering from two or more scattering centers. This is justified 
to leading order in as long as the locations of the scattering centers can 
be considered random on the scale of Af. In other words, the interference can 
be neglected in an ensemble average if for each scattering center j we have a 
distribution function Tij(r) (J d 3 rrij(r) = 1) that is smooth on the scale of Af. 
Note that rij(r) can still be a good approximation to the delta function when 
looked at on the scale of £o- ln the case that all scattering centers arc identical, 
it is sufficient to specify only one distribution n(r), which is normalized to the 
total number N of scattering centers: J d 3 rn(r) — N 

In a small neighborhood of a scattering center, the state of the system can be 
strongly different from the bulk. For example, liquid 3 He may be in a solid-like 
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state because of the van der Waals interaction with an impurity. In order to allow 
such changes, the treatment within a scattering center must be fully quantum 
mechanical. However, such a calculation is not possible in practice because even 
the wave function of a quasiparticle is not known in any detail. Therefore, one 
has to introduce the properties of a scattering center by some phenomenological 
parameters. A sufficient description is to know the scattering T matrix or, in case 
of an isotropic scattering center, the scattering phase shifts 5i. It is important 
to notice that these are parameters that appear in the normal state, and they 
can in principle be measured without need to go into the superfluid state. 

The discussion above was limited to the neighborhood of the Fermi surface, 
where the quasiparticle picture is valid. Because of their strong potential, the 
impurities also have an effect outside of this range. As a consequence the pairing 
interaction, for example, is modified by the impurities. However, this effect is 
of short range because of the relatively high energy. Therefore we can neglect 
such effects relying on the assumption that the volume fraction of the scattering 
centers is small. Similar argument gives that the Landau Fermi-liquid parameters 
(F?, F*), as well as the density (p) and the dipole-dipole interaction constant 
(<jr<i) in 3 He are unchanged. 

The size of a scattering center has to be small compared to the length scale 
one is interested to study. In the superfluid state this scale is typically set by the 
coherence length £ - Thus a scattering center has to be small compared to £ , 
although it can be large compared to Af. It should be noted that the size of the 
scattering center does not limit the size of the impurity. In order to consider a 
macroscopic body, we represent its surface by a set of scattering centers. This is 
possible because for all practical bodies, the height-height correlations (C(r)£(r+ 
8r)) of the surface have a distribution that is much wider that Af when the 
separation Sr is on the order of £ - Thus the scattering becomes incoherent 
on this scale, and therefore can be represented by different scattering centers. 
Assuming that the particles cannot penetrate into the body, the volume of the 
body should be excluded from the calculation. For example, we estimate the 
volume of scattering centers in aerogel as the volume of few atomic layers on 
the strands (w XfA) rather than the volume of the strands themselves (V). (It 
happens, however, that both quantities are on the same order of magnitude in 
this case). 

Let us comment the effect of the pairing state. Superfluidity arises because 
part of the quasiparticles are weakly bound to pairs. The momenta of the quasi- 
particles in a pair can be denoted by p and — p because the total momentum 
of the pair is small. The dependence of the pair wave function on the direction 
P = p/p can be described using spherical harmonic functions and classified as 
s, p, d, etc. states. In most superconducting metals the Cooper pairs form dom- 
inantly in an s-wave state. In 3 He the pairs form in a p-wave state, and there is 
increasing evidence that d waves are dominant in high-T c cuprate superconduc- 
tors. The amplitude of the superfluid state and the transition temperature can 
be calculated for these cases as a function of impurity scattering, as will be dis- 
cussed in more detail later. In the quasiclassical approximation one finds that the 
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impurities have no effect on a homogeneous s-wave superfluid. For other waves 
the superfluidity is suppressed by impurity, and it disappears when the mean 
free path is reduced to ~ £o- The interpretation is that scattering changes the 
momenta of the quasiparticles, but this has no effect in the s-wave case because 
the wave function is independent of the momentum direction. For other waves 
the sign of the wave function changes for certain scattering angles, which leads 
to destructive interference. There is an effect for s waves, too, if the superfluid 
state is inhomogeneous. Here the scattering makes the quasiparticles more local- 
ized, and therefore they see more the same order parameter. This effect makes 
the superconductor less sensitive to magnetic field, i.e. causes the change from 
type I to type II superconductivity, as mentioned in the introduction. 

All the discussion above has been in weak-coupling limit. This means neglect- 
ing all corrections that are proportional to the small ratio Af/£o- This approx- 
imation is clearly inadequate for some properties of 3 He. For example, it gives 
that the B phase is always more stable than the A phase, contrary to the phase 
diagram in Fig. 1. There are some calculations of strong coupling corrections 
in the pure 3 He (Serene and Rainer 1983). It seems very difficult to take these 
corrections into account in the impure case, and therefore we do not consider 
them here. 

4 Homogeneous Scattering Model 

We do not know the precise distribution n(r) of scattering centers in aerogel. 
Therefore we have to take a model for n(r). The simplest possible one is the 
homogeneous scattering model (HSM). There one assumes that n(r) = n is a 
constant. In other words, the probability for a quasiparticle to be scattered is 
the same at all locations r. This same approximation is commonly used to study 
impurities in superconductors (Gorkov 1959, Abrikosov and Gorkov 1961). 

The mean free path in the HSM is given by I = (ncr) -1 . Here a is the scatter- 
ing cross section of a scattering center. We additionally assume that the medium 
is isotropic, i.e. I is independent of the direction of quasiparticle momentum. We 
also assume that the scattering is nonmagnetic. This means that the scatter- 
ing probability is the same for both directions of the spin, and the spin is not 
changed in the scattering. 

A convenient property of the isotropic HSM is that both the Ginzburg- 
Landau (GL) theory and Leggett's theory of NMR (Leggett 1974) have the same 
form as in pure 3 He. The changes appear only via modified parameter values of 
these theories. The GL theory allows a convenient way to represent the results 
of the HSM at temperatures near T c . Therefore we describe it in detail below. 

The order parameter in superfluid 3 He is a complex 3x3 matrix A^. This 
describes the wave function of the Cooper pairs such that for a given direction p 
of momentum, the (unnormalized) spin wave function is given by (Leggett 1975) 

VVn(p) = {-d x + id y )\ TT) +d z \U + IT) + (d x + id y )\ ||) , (2) 
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where = J2i A^Pi. This apparently strange notation has the advantage that 
the spin \i and the orbital i indices in A^ transform equally in coordinate ro- 
tations. The amplitude of A^ describes how strong the superfluid state is. In 
particular, it goes to zero continuously when the temperature approaches the 
superfluid transition temperature T c . 

The free energy / of the superfluid is a function of the order parameter A^. 
Similar to the Ginzburg-Landau theory of superconductivity, one can expand / 
in powers of A^i when the temperature is not much different from T c . The most 
important "bulk" terms in the expansion are (Mermin and Stare 1973) 

/bulk = fn + aA^A^ + PM^A^A 2 + (3 2 (A lll A*^) 2 

+foA fli A vi A V jA ll j -\- PiA^AviAyjApj -\- fisApiAviAvjApj . (3) 

Here a summation over repeated indices is implied. Equation (3) includes the 
allowed terms up to fourth order, since several terms have to vanish by sym- 
metry: / has to be real valued, and it should remain unchanged in rotations 
of both the spin and the orbital parts of A^. This is because the pairing in- 
teraction is unchanged in such rotations. For stability the fourth order terms 
have to be positive definite, but otherwise the coefficients p\ are arbitrary in a 
phenomenological approach. 

Minimizing (3) one finds the normal state (A^ = 0) when a > 0. In the 
superfluid state (a < 0) there are several possible minima. The most important 
are the ones corresponding to the A and B phases. The order parameter in the 
B phase has the form 

A M = cxp(iti)A B R M , (4) 

where R^i is an arbitrary rotation matrix (R^iR^j = i? an arbitrary phase, 
and 

p r 3 2 3a 

?b ~ U = -zaA B = - — -. (5) 

2 4(3/?i2 + P345) 

We use the notation p^... = + (3j + The A phase has 

A-nj = A A d^(rhj + ifij) , (6) 

where d is an arbitrary unit vector and m, n, and 1 = m x n from an arbitrary 
orthonormal triad. The energy and amplitude are 

/a - /„ = aA 2 A = . (7) 

4P245 

In order to get the coefficients a and p\ one has to solve the quasiclassical 
equations. This is briefly explained in the Appendix. The result for a is 



N(0) 
a = — - — 



Ten ' ^— ' V n — O n ~ O 
n—1 v * * 



(8) 



where 2N(0) is the density of states at the Fermi surface and x = hvf /4:irT£ tI . 
The scattering parameter is the transport mean free path £ tI = (notr) -1 - The 
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Fig. 2. The transition temperature in aerogel relative to that in bulk, T c /T c o. The hor- 
izontal axis is the coherence length £o = hvf/2nkBT c o- The experimentally inaccessible 
values of £o are denoted by shading. The data points are the same as in Fig. 1: Porto 
and Parpia (1995, A), Sprague et al. (1995, ■), and Matsumoto et al. (1997, o). The 
solid lines correspond to the homogeneous scattering model (HSM) at mean free paths 
£ tr = 320 nm (a) and 230 nm (b), and the dashed lines correspond to slabs of thick- 
nesses D — 105 nm (c) and 74 nm (d). The dotted line is the isotropic inhomogeneous 
scattering model (IISM) with average mean free path 70 nm, unit cell radius R = 140 
nm, and scattering profile parameter j = 8. 



transport cross section atr gives more weight to large scattering angles 9 that the 
total cross section a. The precise definitions are atr = 47r((l— cos 9) (da / dQ) (9)) n 
and a = 4n((da / 'df2)(6)) q. Here (da / dQ)(9) is the differential scattering cross 
section and (. . .)q denotes average over 4-7T solid angle. 

The a term (8) was first calculated by Larkin (1965). It has exactly the 
same form as calculated earlier for magnetic impurities in s-wave superfluid 
(Abrikosov and Gorkov 1961). The transition temperature T c obtained from the 
condition a(T c ) = is plotted in Fig. 2. The figure also shows the experimental 
data plotted as a function of the coherence length £o- We see that the HSM 
gives a suppression of T c , but its dependence on £ is much weaker than is seen 
experimentally. 

We have assumed in Fig. 2 that £ t r is a pressure- independent constant. We can 
identify three effects that could make £ tI — (na tI )~ 1 pressure dependent, (i) The 
structure of aerogel changes, (ii) The Fermi momentum changes, (iii) The wave 
function of a quasiparticle changes. It seems that all these effects are negligible 
compared to the changes shown in Fig. 2. In particular, the Fermi momentum 
hkf changes only 10% over the whole pressure range, and calculations with hard 
spheres indicate that <7t r is very weakly dependent on k{ for impurities that are 
large compared to Af . What remains most uncertain is point (ii) because we do 
not know the wave function of a quasiparticle and how it changes with pressure. 

In the neighborhood of the transition temperature T c we expand a(T) = 
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a 1 x (T/T c — 1). The higher order corrections in the Taylor series are dropped 
because they would not improve the accuracy of the GL functional (3). From 
equation (8) one easily gets 



a = 



N(0) 



l-x c ^(n- \ +x c )~ 



n=l 



(9) 



where x c — hvf / 4irT c £ tr . 

It is not easy to calculate the fourth order coefficients /3, with the same 
generality as a. Therefore we make the additional assumption that only the 
s-wave scattering phase shift <5o is nonzero. In this case <7t r = <J, and we get 



ft 

W 



a = 



V -i J 

N(0) 





1 

\-lJ 



(Apf\ 
1 Apf 1 

A/% 

\ApfJ 



(10) 



N(0)hv f 2 „^ 



k + x c y 



n=l 



Here Z\/3| c are strong-coupling corrections which are not evaluated here. A model 
calculation for them in pure 3 He is given by Sauls and Serene (1980). 

Note that a and a, instead of depending separately on n and the cross section 
a = (47r/fcf) sin 2 So, depend only on the combination I = (an)^ 1 . This is not 
the case for 6, which is also a function of So- The limiting cases sin 2 So — > and 
sin 2 So — 1 are known as Born and unitarity limits, respectively. This additional 
degree of freedom decreases the precision of the theoretical results because So 
is not generally known. However, we argue in the following that this degree of 
freedom is averaged out in aerogel. Firstly, if the impurities are not identical, we 
can expect that there is a distribution of <5o's instead of a single value. Secondly, 
the calculation including all partial waves has been done in the case of B phase 
(Thuncberg et al. 1981). There one finds a similar spread of the results as a 
function of the phase shifts. But for large impurities the phase shifts Si with 
different I are essentially random, as demonstrated by Si for hard spheres. Both 
these arguments support that a reasonable guess, if no more detailed information 
exists, is to assume random phase shifts sin 2 So — » 0.5. 

A necessary condition for the stability of the A phase is /a < /b- Using 
equations (5), (7), and (10) this reduces to 



a < QA(3{ C + 2A/3f - 4Af3[ 



45 



(11) 



Because a increases monotonically by factor 3 with decreasing T c , the B phase 
becomes more favored with increasing scattering assuming that Z\/3| c do not 
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grow even more. Assuming A(3j c remain constants, we can expect stable A phase 
only at temperatures (in mK) and pressures where it is stable in the pure case. 
Calculations with other phases indicate that they arc not serious competitors to 
A and B phases. 

In order to make more comparisons with experiments, we need to consider 
additional terms in the GL functional F = J d 3 rf(r). These are the gradient 
energy 

f k = K[( 1 -l)d i A fli d j A; j +d i A tlj d i A; j ] , (12) 
the energy of the magnetic field H, 

/, - -\xn 'H 2 + g z H^AUH, , (13) 

where Xn is the susceptibility of the normal state, and the magnetic dipolc-dipolc 
interaction energy 

/d - 9*{\Au | 2 + AijA* }i - lA^) . (14) 

Here /k is purely phenomenological similar to the bulk energy (3), but in f z and 
/d we have dropped some terms, which could be there in a pure phenomenological 
theory. For the coefficients we calculate (Appendix) 



K = 



N(0)K 2 vf " 3 



h 2 j 2 N(0) 



(16) 



487r2Tc2(1 + F a )2 E(--|+^)- 3 (17) 



■ 9d= 40 



<=CUt off/27TT c 

n=l 



2 



(18) 



where 7 is the gyromagnetic ratio, Fq a Fermi-liquid parameter, R a renormal- 
ization constant for the dipole energy, and e cu t off a high-energy cut-off (Lcggctt 
1974). The dipole-dipole coupling constant is different form the other coeffi- 
cients [(8)-(10), (15)-(17)] because it is dominated by high-energy processes (it 
would diverge for e C utoff — » 00), and therefore it is a constant (independent of 
scattering) . 

All the coefficients above reduce to the well know GL coefficients in the limit 
vanishing scattering, £ — > 00 (Fetter 1975, Thuneberg 1987). 

In most cases the terms f z and /<j are small compared to /bulk- Treating them 
as perturbations, we find in the A phase (apart from constants) 

f z A^~^Xn 1 H 2 + g z A 2 A (d-H) 2 (19) 
JdA = -2. 9d Zii(d- 1) 2 , (20) 
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which favor dlH and d || 1, respectively. The B phase is slightly more com- 
plicated. We parameterize the rotation matrix by an angle 9 and an axis n 
of rotation. The dipole-dipole energy /d favors 9 = arccos(— 1/4) = 104°. The 
external field leads to distortion of B phase order parameter (4). This gives rise 
to the energy term 

/ d zB = -g^(n.H) 2 , (21) 

4P345 

which favors n || H. 

For comparison to experiments, we still need to relate the observables to the 
GL coefficients. The torsional oscillator experiments measure the superfluid den- 
sity p s (Porto and Parpia 1995, Matsumoto et al. 1997). This can be calculated 
by evaluating /k (12) in the presence of a phase gradient: Vj4 w - = iqA^j, where 
the superfluid velocity v s = hc\/2mz and 7713 is the mass of a 3 He atom. The main 
quantities measured in the NMR experiment are the magnetic susceptibility \ 
and the shift of the resonance frequency Su> from the Larmor value. The former 
is obtained by evaluating / z . The latter needs evaluating /d in the dynamic case 
(Leggett 1974). The results for homogeneous A and B phases are 



p sA = -^KA 2 A [>y + 1 - (7 - 1)(1 • v s ) 2 ] (22) 
f^ = ^(^2)KAi (28) 



8m 2 

¥ 

8m 2 

r? 

XA = Xn (dl H) (24) 

Xb = Xn - 2g z A 2 B (25) 

5u; A = ^-g d Al [(i • d) 2 - (i • H) 2 ] (dlH) (26) 

5uj b = ^Lg d A 2 B \h x H| 2 (9 = 104°) . (27) 

Here v s and H denote the directions of the superfluid velocity and field, respec- 
tively. The susceptibility of the normal phase Xn has to include the inert layer 
of 3 He atoms on the aerogel strands, which in fact is the dominant contribution 
at low temperatures (Sprague et al. 1995). The frequency shifts 5tu are for small 
tipping angles of the magnetization from the equilibrium direction H, and the 
magnetic field H is assumed large in comparison to y #d/ffz ~ 0.2 mT. The re- 
sults for xa, 5cja, and Su>b are limited to the conditions indicated in parenthesis 
after each equation. 

The NMR experiments with 3 He (no 4 He mixed) see no deviation of the sus- 
ceptibility from Xn (Sprague et al. 1995, 1996). In the HSM model, we interpret 
this as evidence for the A phase. Therefore, the measured frequency shift is best 
compared with Su>a- In Fig. 3(a) we plot the suppression factor 

= XnM*Tc) _ A 2 (tT c ) _ 
X " SUA ~ XnO<Wo(tT c0 ) - A 2 (tT c0 ) ~ bA2 ■ (2 * } 
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Fig. 3. The suppression factors for A 2 and p s as a function of T c reduction: (T c /T c o)' 2 . 
The lines are theoretical results for the B phase when t — > 1 [see equations (28) and 
(29)]. The dashed lines denote Born and unitarity limits, and solid lines are for the 
intermediate case sin 2 So = 0.5. The data points are from Sprague et al. (1995, A type 
phase, ■), Sprague et al. (1996, B type phase, □), and Porto and Parpia (1995, A), all 
corresponding to t « 0.9. 



Generally, the suppression factor is the ratio of a quantity in aerogel relative 
to the same quantity in the pure case. The latter is denoted by subindex 0. The 
ratio is taken at the same temperature t measured relative to the corresponding 
T c . We consider here t in the range 0.9 ... 1 where the GL theory is expected to 
be valid. The middle equality (28) is based on equation (26). It is valid because 
<7d is a constant and the minimum-energy orientations of d and 1 are unchanged 
by scattering. (The latter is not necessarily satisfied in generalizations of the 
HSM, see Section 6.) Similarly we can define a suppression factor for p s : 



The B phase p s o is used as a reference in the experimental points of Fig. 3(b). 

The theoretical suppression factors S^2 and S Ps are plotted in Fig. 3. In these 
the impure B phase is compared with bulk B phase. Very similar suppression 
factors are obtained when impure A phase is compared with bulk A phase: The 
solid lines are the same for A 2 and also for p s (22) averaged over all orientations 
of the A phase [corresponding to (1 • v s ) 2 — ► |]. The variation between Born and 
unitarity limits is slightly smaller in the A than in the B phase. 

We see from Fig. 3 that HSM indeed gives a suppression of A 2 and p s , but 
it is insufficient to explain quantitatively the measured suppression. 




(29) 
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5 Inhomogeneous scattering models 

The disagreement between the experiments and the HSM is clear in Fig. 3. We 
have checked that this failure does not arise from limitation to the neighborhood 
of T c or to the s-wave scattering approximation. Also magnetic scattering does 
not seem to give a solution to this problem. We also discuss in Section 6 that 
the problem cannot be explained away by random anisotropy of the aerogel. 
Therefore it seems inevitable to sacrifice to basic assumption of the HSM: the 
homogeneity of the scattering. 

A simple model of inhomogeneous scattering is to consider 3 He between two 
diffusively scattering planes. The transition temperature for this "slab model" is 
calculated by Kjaldman et al. (1978), and the suppression factors are evaluated 
by Thuneberg et al. (1996). This gave promising results (see Fig. 2) but it is 
not a suitable model for aerogel because of its strong anisotropy. Recently an 
"isotropic inhomogeneous scattering model" (IISM) was studied (Thuneberg et 
al. 1997). This is consistent with the observed isotropy of aerogel, and it gives 
rather good fits to both T c and Because only preliminary calculations on 
the IISM has been done, we leave its discussion to another occasion. What is 
of importance here that part of the HSM seems to remain valid: Although the 
gap amplitude A 2 is badly overestimated in the HSM, the equations (22)-(27) 
for p s , x, an d 5lo may still constitute a reasonable approximation. This will be 
used in the following section. Similar arguments may be applied to the HSM 
calculations by Baramidze et al. (1996). 



6 Anisotropic HSM 

In this section we consider how to generalize the isotropic HSM (Section 4) 
to anisotropic scattering. The anisotropy will affect all the terms in the GL 
functional, but the most important effect comes from the modification of the 
second order bulk term (3). The anisotropy requires the replacement 



n,A* A ■ 



ni- ■ A* A ■ 



(30) 



So the scalar a is replaced by a tensor a. Solving the quasiclassical equations in 
this case gives 
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Unit matrices multiplying m — ^, for example, are not shown explicitly. The 
transport cross-section tensor is defined by the double angular average 
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(32) 
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where da/dfi is the differential scattering cross section and the total cross section 
<r(k) = 4ir((d<j/df2)(lg., k')}^,. [Inequation (31) we have assumed that both k and 
k' dependencies of (da/df2)(k, k') can be represented by s and p-wave spherical 
harmonic functions.] For estimation of the anisotropy in aerogel we consider a 
long rod, which scatters diffusely (randomly) the particles hitting it. This gives 
na tT = £ _1 (| — |aa), where a is the direction of the rod and I denotes the average 
(transport) mean free path. Assuming the anisotropy is small, the leading effect 
of the anisotropy can be represented by 

/a = -ax&iA^A^aj , (33) 

where 

<*-^f>- i+ ..)-. (34) 

m— 1 

It is crucial that the anisotropy direction a is not a constant. We expect that the 
orientational correlation decays in a length L a ~ 20 nm, which is on the order 
of the average distance between strands. 

The anisotropy (33) shifts the transition temperature. It also changes the 
relative stability of different phases. Instead of discussing these, we will here 
concentrate on how the anisotropy can crucially modify the NMR properties. 
Some of the ideas below are suggested by Volovik (1996). 

Quite generally, we can write the following hydrodynamic free energy for the 
A phase 

/a = -iA d (d • l) 2 + ±A z (d • H) 2 + ±A a (a • i) 2 + ± P± v s 2 + ±( P || - Pl _)(\ ■ v s ) 2 

+cv s ■ v x i - c (i ■ v B )(i ■ v x i) + \k s (v ■ Xf + \K t (\ • v x i) 2 
+ix b |ix (vxi)| 2 + ix 5 |(i.v)d| 2 + i^ 6 [(ix v),d,)] 2 . (35) 

The first two terms are already familiar from equations (19) and (20). The third 
is the anisotropy term (33) with A a = 2a\A\. The rest arises from the gra- 
dient energy (12) when the A-phase order parameter (6) is substituted into it 
(Cross 1975). All the coefficients are proportional to A\. Because this common 
factor drops out in relative comparisons of the terms, our estimations below are 
independent of A\. 

The idea is that the random field a tries to orient the 1 vector. However, the 
gradient energy strongly limits the variation of 1 on the scale L a . Instead, 1 varies 
only on a "orbital" scale L Q ^> L a . In addition, the d vector varies on a "spin" 
scale L s , which also is large in comparison to "aerogel" scale L a . We estimate L 
and L s using arguments presented by Imry and Ma (1975). The exact functional 
(35) is approximated by 
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This is a functional of only two variables, L s and L a . The first two terms consti- 
tute an interpolation of the dipole-dipole energy between two limits. The min- 
imum energy — ^Ad is obtained in the limit L s = L Q . In the opposite limit 
L s 3> L Q , the dipole-dipole energy consists of a random average — gAd plus 
a fluctuation —^rj\d(L /L s ) 3 ^ 2 (Imry and Ma 1975). The origin of the factor 
(L /L s ) 3 / 2 is that in a box of volume L z there are N = (L /L s ) 3 uncorrelated 
regions of 1, and therefore the total fluctuation they produce is proportional \/N. 
Here rj = 2/3\/5 is the standard deviation of (n • z) 2 for a random unit vector n. 
The third and fourth terms in (36) approximate similarly the anisotropy energy 
^A a (a • l) 2 , but there only the limiting form L Q >• L a is needed. The last two 

terms are simplified gradient energies of d and 1. We assume K s s=s \2KA\ and 
K Q « 6KA\, where a factor 3 arises from the dimensionality of the space. 

The simplified energy functional can now be trivially minimized numerically, 
and also to a great extent analytically. We find two basically different solutions, 
which are described below. These solutions have equal energy if L a = 48 nm, 
which is only slightly larger than the estimated distance between stands 20 nm. 
The numerical estimates use this and the pressure of 2.8 MPa. 

In the "dipole locked" solution d and 1 follow closely each other, d(r) s=a l(r). 
We find rather long L s m L Q w 70 /im. The NMR experiments are done in large 
magnetic field, so that d w 1 is allowed to vary only in the plane perpendicular 
to H. We sec from equation (26) that the NMR frequency shift Scoa is essentially 
unchanged from the HSM prediction. This state has lowest energy for L a < 48 
nm. 

In the "dipole unlocked" solution 1 varies on a much shorter scale than d, 
L Q <C L s . We find L Q w 8 /im and L s w 3 mm. Again the magnetic field limits d 
to the plane perpendicular to H. However, 1 is free to vary in all directions. The 
NMR frequency shift is totally suppressed: it is reduced by factor 10 -5 relative 
to the HSM. This is because ((1 • d) 2 ) w ((1 • H) 2 ) w | in the expression for 5lu a 
(26). This state has lowest energy if L a > 48 nm. 

As discussed above, the experiments see that the NMR frequency shift is 
smaller than the HSM prediction. However, they still are on the same order of 
magnitude. This strongly supports the dipole-locked state because in the un- 
locked state the suppression would be much more severe. Moreover, because the 
suppression in the locked state is unchanged from HSM, we reach the conclu- 
sion (stated in Section 5) that random anisotropy cannot explain the difference 
between measurements and the HSM. 

The presence of the dipole-locked state of 3 He in aerogel implies that the 
anisotropy of aerogel is very small: it is equivalent to rods whose directions are 
correlated over distance L a < 48 nm. This invalidates the slab model, because 
the thickness D of the slab is larger that this (see Fig. 2). 

Above we have considered the NMR properties only in the case that the "tip- 
ping angle" (3 between the field and the magnetization is small. When the angle 
is increased, the frequency shift 5u) gradually diminishes. This is similar to bulk 
3 He-A, and it is well understood theoretically (Brinkman and Smith 1975, Fomin 
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1978). However, the frequency shift suddenly disappears at about (3 = 40°, and 
remains zero for all larger angles (Sprague et al. 1995). As suggested by Volovik 
(1996), this striking observation can be associated with the locked-^unlocked 
transition discussed above. At the time of writing this, this is still under study, 
so we postpone the discussion to another occasion. 

7 Conclusion 

We have extensively studied the homogeneous scattering model in the Ginzburg- 
Landau approximation. Although it fails to produce the correct gap amplitude, 
it may successfully applied in the hydrodynamic region. As an example we con- 
sidered NMR in the A phase. In the future it could be used, e.g., for studying 
vortices of 3 He in aerogel. 
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Appendix 

In this appendix we briefly explain how the Ginzburg-Landau coefficients (a, 
0i, etc.) are calculated in the quasiclassical theory. An intermediate quantity in 
the calculation is the quasiclassical 4x4 matrix Green's function <?(k, r, e m ). 
The arguments are the direction of the momentum k, the location r, and the 
Matsubara frequencies e m — TrkBT(2m+ 1). The Green's function is determined 
from the Eilenberger equations 

[ie m T3 - v- p- A,g] + iv F k ■ V r g = (37) 
99 = • (38) 

Here fj denote the Pauli matrices in Nambu space, and [A, B] = AB — BA. For 
more details the reader is referred to the review article by Serene and Rainer 
(1983). The self-consistency equations for the diagonal v and off-diagonal A self 
energies are given in formulas (5.10) of Serene and Rainer (1983). The impurity 
self-energy p(k, r, e m ) equals n(r)i(k, k, r, e TO ), where n(r) is the concentration 
of the scattering centers. The {-matrix of a single scattering center is determined 
by 

f(k, k', r, e m ) = «(k, k') + N(0)(v(k, k")g(k", r, e m )t(k" , k' , r, c m )) t „ . (39) 

Also we need the free energy functional which is given by equation (5.11) of 
Serene and Rainer (1983). 
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The coefficients are calculated by solving Green's function perturbatively 

g = go + g\ + §2 + Here g is the normal state Green's function, and the 

subindex denotes order in A Then one collects terms of order 0, 1, 2, and 3 
in the Eileberger, self-consistency and i-matrix equations, and solves them in 
each order. Substitution to energy functional gives the results a (8) and (10). 
For other coefficients it is sufficient to limit to first order in A, but one has 
to do additional expansions in gradients and in magnetic field. Substitution to 
energy functional gives K (15), 7 (16) and g z (17). The dipole coefficient (18) 
is obtained using the dipole-dipole Hamiltonian (Leggett 1974) in the general 
energy functional (5.6) of Serene and Rainer (1983). 
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